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Using a quantum map version of one-dimensional Anderson model, the localization-delocalization 
transition of quantum diffusion induced by coherent dynamical perturbation is investigated in com¬ 
parison with quantum standard map. Existence of critical phenomena, which depends on the number 
of frequency component M, is demonstrated. Diffusion exponents agree with theoretical prediction 
for the transition, but the critical exponent of the localization length deviates from it with increase 
in the M. The critical power Cc of the normalized perturbation at the transition point remarkably 
decreases as Cc ~ (M — 1)~^. 
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I. INTRODUCTION AND MODELS 

The localization phenomena are persistent and ro¬ 
bust in one-dimensional disordered systems (IDDS) [H- 
[^. It still remains even in two-dimensional disordered 
systems. However, if the dimension d of the disordered 
system is more than 2, the localization becomes unsta¬ 
ble, and the localization-delocalization transition (LDT) 
takes place, and finally an irreversible diffusion sets in 
when we consider quantum diffusion of an initially lo¬ 
calized wavepacket. The critical phenomena of LDT 
have been extensively studied on the basis of the one- 
parameter scaliM theory (OPST) of the localization by 
many authors The recent reviews and develop¬ 

ments for LDT have been given in a commemorative book 
0 , and references therein. 

On the other hand, similar localization phenomena 
were discovered for the quantum kicked rotors (KR) typ¬ 
ically exemplified by the quantum standard map (SM), 
and it can be interpreted as the localization phenomenon 
of a class of IDDS in terms of Maryland transformation 
pH - fTsjl . In this context, the additional dimensionality 
(d — 1) corresponds to the number of the dynamical de¬ 
grees of freedom M applied to the KR, and the LDT in 
IDDS corresponds to the ergodic transition when we con¬ 
sider the dynamically perturbed standard maps. Based 
upon this correspondence, the critical phenomenon of the 
LDT was observed for Cesium atoms in optical lattice 
settings 0 El- 

On the analogy of the standard map’s case, we can ex¬ 
pect that the localization in the IDDS is unstable against 
the dynamical perturbations. We have proposed the de- 
localization scenario that the dynamical perturbation to 
IDDS in general enhances the localization length and re¬ 
store the diffusive motion in a strong perturbation regime 
01 . Considering that electrons are interacting with lat¬ 
tice vibrations, the effect of dynamical perturbation by 
phonon modes is essential. It models the fundamental 
dynamical and deterministic process of the quantum 
electronic motion turning into diffusive one which allows 


time-irreversible kinetic description. 

Here, a basic question arises. Whether the LDT hap¬ 
pens in the IDDS under the interaction with dynamical 
degrees of freedom, and if it happens, how the nature of 
LDT changes with increase in the mode number M. To 
answer this we introduce the ID Anderson map (AM) 
described by the unitary time-evolution operator 

= Q-il^T{p)/2h^-iAU(tm)V{q)/h)^-iAT{p)/2h^ 

for wave function defined on the discrete lattice site 
q{= n), where T{p) = 2cos{p/h) = -I- and 

V{q){= V(n)) is random on-site pot ential uniformly dis¬ 
tributed over the range [—IF, W] The dynamical 

perturbation is modeled by the sinusoidal periodic per¬ 
turbation superposed onto the on-site energy as 

M 

f(f) = ^ +-f^^cos{u}kt), ( 2 ) 

where M and e are the number of the frequency compo¬ 
nent and the strength of the perturbation, respectively. 
The time evolution by the operator Um approximates the 
unitary evolution of the dynamically perturbed Anderson 
model, 

^^ du{n, t) _ ^ i) _|_ f{t)V{n)u{n, t),(3) 

at 

for a short time interval A up to the correction of O(A^), 
and the unperturbed Anderson map (e = 0) retains the 
localization properties of the Anderson model 0 , [^ . 
In addition, the numerical verification of the presence 
of LDT is very hard for the perturbed Anderson model, 
therefore we examine here the perturbed Anderson map 
to explore the presence of LDT, where we take A = 1, 
typically. Note that the strength of the perturbation is 
divided by a/M so as to make the total power of long-time 
average independent of M, and the frequencies are taken 
as incommensurate numbers of 0(1) |2l| . Replacing by 
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T{p) = p^/2 and V{q) = Kcos{q), Eq.(IT|) turns into the 
SM, which exhibits the LDT in the momentum space [l^ . 

In the perturbed AM, both localation and delocaliza¬ 
tion have been observed [Il,[2|. However, the nature 
of the transition from the former to the latter was not 
known, in particular, the presence of critical phenomena 
in the transition process is still unclear. In this paper, we 
numerically investigate the critical nature of the LDT in 
the AM in comparison with the LDT in the SM, which 
can be analyzed by the OPST In particular, we 

are interested in the mode number M dependence of the 
transition, regarding the M with the additional dimen¬ 
sion (d — I) according to the interpretation in the case 
of SM . We increase the effective dimensionality 

{M+1) = d far beyond 3. It will provide a crucial test for 
the mean-field theory of the Anderson transition, which 
regards d = 3 as the lower bound above which the critical 
exponents lose the d—denendency. 


II. SUBDIFFUSIVE PROPERTIES OF 
WAVEPACKET DYNAMICS AT THE 
LOCALIZATION-DELOCALIZATION 
TRANSITION 

For M = 1, the localization length increases exponen¬ 
tially with the coupling strengeth e, but we could not 
confirm the presence of delocalized state. However, if 
M > 2, we can confirm the presence of critical state, 
which evidently borders the delocalizing behavior and 
the localizing one. This will be described closely in the 
present section. 


curves spreading like a fan into the delocalization regime 
(e > Cc) increasing up to the normal diffusion m 2 {t) —>■ 
Dt (A —>■ and localization regime (e < Cc) decreas¬ 

ing down to the localization length m 2 {t) —^ ^ (A —t““). 
Here, the diffusion constant D = D{e) and the localiza¬ 
tion length ^ = ^(e) approach to zero and infinity, re¬ 
spectively, as e —Cc- 



A. Numerical results 

Let us introduce the on-site probability P{n,t) = 
\u{n, t)p. The main tool of our analysis we use the time- 
dependent mean square displacement (MSD) m 2 {t) =< 
>YP{n-,t) >n of the propagating 
wavepacket starting from localized one, u{n,t = 0) = 
Sn,noj where < • • • >n denotes the ensemble average over 
different random configuration of V (n). With increase 
in the perturbation strength, the time-evolution changes 
from the localized behavior to delocalized one passing 
through the critical behavior at a certain critical strength 
e = Cc. (See Figllja) for the result of AM with M = 2. 
) At the critical state the MSD exhibits a power-law 
asymptotic (subdiffusive) dependence m 2 {t) ~ char¬ 
acterized by the diffusion exponent a, which is close to 
the theoretical value 0.66 discussed later. Instead of the 
MSD, we introduce the scaled MSD 

( 4 ) 

with respect to the critical behavior and show its tem¬ 
poral evolution at various e including Cc in Figlljb). The 
critical curve indicated by the bold line separates the 


FIG. 1: (Color online) The diffusive properties of the 
wavepacket in the perturbed AM with M = 2. The system 
and ensemble sizes are N = 2^® ~ 2^^ and 10 ~ 100, respec¬ 
tively, throughout this paper, and we mainly take W — 0.5 
as the disorder strength and ft = 0.125 as the Planck con¬ 
stant, respectively. (a)The double logarithmic plot of m 2 (t) 
as a function of time for different values of the perturbation 
strength e, where the diffusion exponent a(— 0.66) is deter¬ 
mined by the least-square-fit for m 2 {t) with the critical case. 
(b)The scaled MSD A(e, t) as a function of time for different 
values of the perturbation strength e. Note that this is log- 
log plot. {c)Ngp as a function of time for different values of 
the perturbation strength t. Note that the black thick lines 
in the panels (a),(b) and (c) show the results at the critical 
case Cc = 0.075. The blue dashed curves show the results for 
e < Ec in the panels (a), (b) and (c). (d)Semi-log plots of 
the scaled probability density P{ns{i)) as a function of the 
Us = n/ ^ 7712 ( 1 ) for t = 4 X 10^, 12 x 10^, 20 x 10^ at the 
critical case. The curves of all cases are well-overlapped. 


To clarify the shape of the distribution, in addition 
to the MSD, we introduce the non-Gaussian parameter 
(NGP) Nqp defined by 


Nap{t) 


1 7714 (t) 

3 m2{tY 


( 5 ) 
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where m 4 {t) = ^{n— < n >)‘^P{n,t). Figure [TJc) de¬ 
picts the time dependence of NGP for various e. At the 
critical point the NGP keeps the same nonzero-value, im¬ 
plying that the shape of the distribution function takes 
a similar non-Gaussian form throughout the time evolu¬ 
tion [H [ 23 . Figure [ijd) shows the distribution func¬ 
tion Ps{ns{t),t) = P{n,t)dn/dns{t) at several Ps as a 
function of scaled coordinate ns{t) by the spread of the 
wavepacket for AM with M = 2 as 


ns{t) 


n n 

- ; OC -777. 


( 6 ) 


B. Comparison with theoretical prediction 

According to the mean-field theory of the Anderson 
transition in the d—dimensional disordered systems [^ . 
the subdiffusion m 2 {t) ^ appears only at the critical 
point Cc and the exponent is represented by the formula 
ctM = I = for 2 < M < 10. On the other hand, 
the exponent (3 has been supposed to be related with a. 
For example, a phenomenological theory based upon the 
assumption that the distribution function is described 
by a master equation with memory kernel tells that the 
relation between a and /3 is an universal relation 


Evidently, the scaled representation Ps{ns,t) does not 
have explicit t dependence as is expected. Thus we de¬ 
note the scaled distribution function simply by Ps{ns{t)). 




FIG. 2: (Color online) (a)The scaled invariant distribution 
P{ns) of AM as the function of ris = n/-^ 7712 ( 1 ) for M — 
2, 3, 4 at each critical perturbation strength Cc from the inside 
to the outside. (b)The plot of In | — lnP(ns)| as a function 
of |ns| in the logarithmic scale. The slopes correspond to the 
exponent /3 of the stretched Gaussian distribution. 


We further investigate the invariant function form of 
the wavepacket at each critical point of various M’s, as 
seen in Figl2](a) for AM with M = 2, 3,4. It is suggested 
that the tail of the scaled invariant shape of the distribu¬ 
tion function takes the stretched Gaussian distribution 

P{ns)exp{-\ns{t)f), (7) 

except for the range close to the origin of the critical 
state, where P is the distribution exponent (streched 
Gaussian exponent). The tails are shown in[2](b) in the 
plot of In I — In P(ns) I as a function of In | rig | for each case 
in FigHKa). The slopes of the plots correspond to the ex¬ 
ponents /3 of the stretched Gaussian distribution, which 
are decided by the least-square-fit except for the range 
close to the origin. The M-dependence of estimated dif¬ 
fusion exponent a and the stretched Gaussian exponent 
j3 are summarized in Figl^l 


for 0<a<l,l<,5<2 [23|. This relation predicts 
Pm = {M + 1)/M if the mean-field result is applied. The 
results of our numerical experiment are compared with 
the theoretical prediction in FiglSl and the data almost 
agree with the relation over a wide range 2 < M < 10. 
(The inset of FigUlshows the a—P relation.) We executed 
the same analysis for SM, concluding the results for AM 
has a great deal in common with those of SM for 2 < 
M < 10, as seen in FiglH] 



FIG. 3: (Golor online) The diffusion exponent a and distri¬ 
bution exponent /? as a function of M for AM and SM. The 
broken lines are theoretical predictions om and Pm, respec¬ 
tively. The inset shows the plot of /I as a function of a. The 
points {a = 0, P = 1) and (a = 2/3, P = 3/2) are denoted by 
filled circles which correspond to the exponential localization 
and subdiffusion at d = 3, respectively. The green solid curve 
is the theoretically expected universal relation. The parame¬ 
ters are same to the case of Fig[T]for AM, and N = 2^“^ ~ 2 ^®, 
K = 3.1, h = 222^ are for SM. 

It is found that the numerical results almost agree with 
the theoretical predictions except for the stretched Gaus¬ 
sian exponent P in the perturbed SM with the larger 
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effective dimension (M + 1) > 8. It seems that the dis¬ 
agreement of P between the theoretical prediction and 
numerical data is due to the insufficiency of the ensemble 
and system size for the tail of the invariant distribution 
functions. 

III. FINITE-TIME SCALING ANALYSIS OF 
THE LOCALIZATION-DELOCALIZATION 
TRANSITIONS 

We next investigate the critical exponent v related to 
the localization (correlation) length which is supposed 
to diverge ^ ~ |e — for the localized regime e < Cc 
(for the diffusive regime e > Cc). The LDT can be in 
general observed both for AM and SM. We exhibit here 
the finite-time scaling analysis for AM taking the case of 
M = 5 as the example. First, we show in Fig0Kb) the 
observed value of lnA(e,t) at various different times tm 
as a function of e. A remarkable feature is that all the 
curves crosses at a single point, which can be regarded 
as Cc- This fact allows us to follow the OPST which is 
usually supposed for the LDT as follows: 

A(e,t) = F^((e, - ( 9 ) 

where F{x) is a differentiable scaling function. We 
note that the asymptotic function form of F{x) should 
be F{x) —>■ \x\~‘^'^ in order to represent the localiza¬ 
tion A(e,t) —>• t°‘x{e) ~ t °‘\e — With the 

above hypothesis, the scaled MSD can be expressed by 
lnA(e,f) — InAc(t) oc (cc — around the critical 

point e = Cc, where Ac = F'(O). Applying this relation 
to the curves in FigdKb), we can plot tm and the cor¬ 
responding slope s(tm) = ln(A(e,tm)/Ac(tm))/(ec - e) , 
which should oc as shown in the inset of Fig0Ka), 

which enables to decide the unknown exponent v by us¬ 
ing the known diffusion exponent a. With this v, we 
can explicitly construct the localization length function 
^s(e) = ■Co(ec — and further the scaling function 
F{^s{e)t°‘^’^) as a function of a; = ■ Figure 

SKa) shows the scaling functions F{x) constructed by the 
time-dependent data of MSD at various different e’s close 
to Cc- All the time-dependent data obtained at different 
e form a unified function if the scaled variable x is used, 
which proves numerically the validity of the OPST. 

In FigHKc), we compare the localization length func¬ 
tion ^s(e) decided indirectly by OPST in the critical re¬ 
gion e ^ Cc with ^n(e) decided directly by the saturated 
MSD data which are precisely calculatable for e’s much 
less than the critical region. The e—dependence of these 
two localization lengths, ^s(e), seem to connect 

continuously, which implies unexpected wideness of the 
critical region in which the OPST works. Accordingly, 
Eq.® based on the OPST immediately leads to 

m 2 {t) oc (e — ec)~^''x^''F{x), (10) 

where x = {tc — and F[x) —>■ x~'^'' for a; —>■ 00 . 

It describes the universal relaxation process toward the 



In X 




FIG. 4: (Color online) The results of the critical scaling anal¬ 
ysis for the perturbed AM. (a)The scaled MSD A(e, t) as a 
function of a; = for some values of e in the AM of 

M = 5. The e-dependent localization length 5s(e) is deter¬ 
mined by a scaling relation Eq. by least-square fit in the 
inset of the panel (a). (b)The same data as panel (a) but plot¬ 
ted as a function of the perturbation strength e. In the ideal 
case, all lines have a common crossing point at e^ = 0.0175. 
(c)The localization length ^(e) as a function of (cc — e) for 
M — 2,3, 4. The filled symbols denote the numerical data di¬ 
rectly obtained by ^m 2 {t —>■ 00 ) in the long-time limit. The 
open symbols indicate the localization length 5s(e) obtained 
by OPST in the critical region. Note that the axes are in 
logarithmic scale. 


localized state, but we do not still know to what extent 
the universal relaxation dynamics works out of the criti¬ 
cal region. 


IV. CRITICAL EXPONENT AND 
PERTURBATION STRENGTH OF THE 
LOCALIZATION-DELOCALIZATION 
TRANSITIONS 

In FiglSKa), we compare the results of iz for AM and 
SM at various M in comparison with theoretical predic- 
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FIG. 5: (Color online) (a)The dimensionality (M + 1) = d 
dependence of the critical exponent v which characterizes the 
critical dynamics in the AM and SM. The red solid line and 
green dashed line are the results of the analytical prediction by 
vvw and vq, respectively. Thick line denotes the lower bound 
by the Harris’ critical inequality. (b)The critical perturbation 
strength as a function of (M — 1) for the AM and SM. Note 
that the axes are in logarithmic scale. Here, we used the SM 
with K = 3.1, ft= 222^11. 


tions. The critical exponent v obtained from the self- 
consistent mean-held theory of the localization (VW- 
theory) is ww = l/((i—2) for 2 < d < 4 and vyw = 1/2 
for d > 4 . The semiclassical theory of Garcia predicts 

1/(5 = 1/2 -|- l/(d — 2 ) which asymptotically approaches 
the value 1/2 of vvw for d —>■ oo m On the other hand, 
the inequality i/ > 2 /d is proposed at the critical point 
by Harris [^. Our results tell that for the larger value 
of M{> 7) the critical exponents of the AM and SM be¬ 
come signihcantly lower than the theoretical lowest value 
1/2 but satisfy the Harris’ inequality [26l - [^ . 

Finally, we show the M-dependence of the critical 
strength Cc in FiglSKb). Note that in the definition of 
e we normalized the perturbation by '/M in order to 
make the power strength of perturbation is independent 
of M for the fixed e. In spite of such a normalization, the 
critical perturbation strength depends strongly upon M. 
As shown in FigEKb), our data for AM and SM indicate 
the inverse power law 

eeCx(M-l)-^ (11) 


up to M = 10. The powers are estimated as 5 ~ 1.1, 
for the AM {W = 0.5 and 0.8), and 5 ~ 1.0 for the SM 
[K = 3.1). Ea. dTT]) means that the total power of the 
perturbation, which is given by is asymptotically 
proportional to ~ 1/M. Such a strong Ad—dependence is 
highly nontrivial and the theoretical derivation has not 
been given to the best of our knowledge. It manifests 
that the LDT with large M is a cooperative phenomenon 
among the degrees of freedom of perturbation and the 
driven system. 

It is quite interesting that the numerical data of Garcia 
and Guevas reporting the delocalization potential thresh¬ 
old of a high-dimensional disordered tight-binding model 
suggests (d— 2 )“^, which seems to be closely related with 
our results 0 . 


V. CONCLUSION 

We investigated critical phenomena of LDT exhibited 
by polychromatically perturbed AM, which models IDDS 
perturbed by coherent dynamical perturbations, in com¬ 
parison with the SM under the same perturbations. We 
confirmed the presence of critical phenomenon for the 
mode number M > 2. The diffusion exponent a and 
distribution exponent (3 agree well with the theoretical 
prediction for M < 10. On the other hand, the criti¬ 
cal exponent v is significantly lower than the predictions 
of the mean-field theory for large M, but it does not 
violate the critical inequality. The critical value of nor¬ 
malized perturbation strength exhibits a remarkable M- 
dependence as Cc (M— 1)“^. Asa result, all the critical 
characteristics of AM agree surprisingly well with SM in 
spite of their fundamental difference. Our results open a 
new possibility of controlling electronic localization and 
conduction by means of externally applied stimulus im¬ 
plemented by optical and/or acoustic devices (^ . 
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